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Two  methods  of  treating  the  buried  insulated  antenna  are  described,  transmission 
line  and  moment  methods.  These  two  methods  are  compared.  Calculations  are 
performed  and  results  given  at  HF  for  a  range  of  antenna  lengths  and  insulation 
thicknesses,  giving  the  gain,  front  to  back  ratio  and  impedance. 
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1  INTRODUCTION 

The  type  of  antenna  considered  is  shown  in  figure  1  and  consists  of  a 
buried  central  conductor  probably  covered  in  a  layer  of  dielectric  insulation 
for  some  part  of  its  length.  This  type  of  antenna  is  described  in 
reference  I  1] .  The  current  can  be  thought  of  as  radiating  an  electric  field 
which  passes  to  the  surface  and  is  refracted  towards  the  interface  on  passing 
through  the  interface.  Normally  most  of  interest  is  the  field  on  the  earths 
surface  and  therefore  also  the  lateral  wave  which  can  be  thought  of  as 
travelling  along  the  interface.  The  angle  at  which  the  rays  travel  to  the 
surface  will  in  this  case  be  the  critical  angle.  Reference  [  2]  suggest  chat 
tilting  the  antenna  at  right  angles  to  these  critical  rays  would  increase  the 
signal.  Due  to  the  difficulty  of  burial  of  long  wires,  only  shallow  burial 
and  antennas  parallel  to  the  interface  are  considered. 

The  central  conductor,  the  insulation  and  the  higher  permittivity  lossy 
ground  form  a  transmission  line  which  will  support  a  travelling  wave. 

Stratton  [  3]  page  547  gives  an  equation  for  the  permissible  values  of  the 
wavenumber  of  the  line  Kl*  If  Che  radius  of  the  insulation  is  much  less  than 
the  wavelength,  and  the  complex  dielectric  constant  changes  appreciably  across 
the  insulator  ground  boundary,  only  the  lowest  order  mode  is  non-evanescant . 
Lee  et  al  [ 1]  give  approximate  formulae  found  in  section  2  for  the  wavenumber 

V 

So  that  there  shall  be  progressive  reinforcement  of  the  wave  at  the 
surface,  the  nearer  is  to  the  wavenumber  K  of  the  lateral  wave  (same  as 
the  wavenumber  in  free  space)  the  better.  The  effect  of  insulation  thickness 
is  dealt  with  in  section  2. 
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Another  factor  of  importance  is  the  front  to  back  ratio.  Although  in 
figure  1  the  waves  are  shown  travelling  to  the  right,  depending  on  the 
character  of  the  current  on  the  central  conductor  there  can  be  a  similar  signal 
to  the  left.  If  the  line  wavenumber  is  fairly  close  to  freespace  then  a  wave 
travelling  from  left  to  right  down  the  line  sets  up  the  phase  of  the  current 
on  the  line  to  favour  a  lateral  wave  travelling  from  left  to  right  ie  the 
phase  of  rays  reaching  the  surface  will  be  that  of  a  wave  travelling  from  left 
to  right.  If  at  the  end  of  the  line  the  travelling  wave  is  reflected  backwards, 
there  will  be  an  appreciable  signal  in  the  backward  direction.  Also  if  the  line 
is  appreciably  less  than  a  wavelength  long  the  backward  signal  will  be  high. 

Another  method  available  for  the  treatment  of  the  problem  will  be  discussed 
in  section  three,  the  moment  method.  In  this  technique  the  buried  antenna  is 
divided  into  elements  and  the  reaction  of  one  element  on  another  in  the  presence 
of  the  interface  above  and  the  dielectric  insulator  is  calculated.  This  method 
allows  a  more  realistic  antenna  configuration  to  include  earth  spikes,  loading 
and  tapering  of  the  insulation  of  the  antenna  as  well  as  including  the  effect 
of  the  interface. 

These  two  methods  allow  design  parameters  such  as  insulation  thickness, 
depth  of  burial  to  be  determined. 

2  TRANSMISSION  LINE  TREATMENT 

As  stated  in  the  introduction  the  central  conductor,  the  insulation  and 
the  outer  conducting  dielectric  (the  earth)  form  a  transmission  line.  The 
equation  on  page  547  of  Strattons  book  [  3]  will  give  the  value  of  the 
wavenumber.  Lee  et  al  [  3]  from  a  different  approach  give  two  approximate 
formulae  for  the  wavenumber  K^:- 


-  K2(l  +  HQ (K^b) / (K^b . ln(b/a)H^ (K^b) ) ) 


1/2 


Accurate  for 


K4/K2 

K2b 


greater  than  16 
much  less  than  1 


\  *  K2((K^(HQ(K4b)  +  K^b.lnCb/ajHjO^b))) 

/(K^O^b)  +  K^b.lnCb/aJH^b)))172. 


Accurate  for 


K4/K2 

K2b 


as  low  as  1.5 
much  less  than  1 


where  K2  is  the  insulation  wavenumber 
is  the  ground  wavenumber 
a  is  the  central  conductor  radius 
b  is  the  insulation  radius 

Hq,  Hj  are  Hankel  functions  of  the  first  kind  of  zero  and  first 
order  respectively. 
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Both  these  formulae  may  be  deduced  from  the  formula  in  Stratton's  book  by 
suitable  approximations.  For  example  by  replacing  the  Bessel  functions 
J<3»  Yo*  Jl»  Y1  by  the  first  terms  in  their  expansions  for  small  arguments* 
formula  1  results.  By  computing  exact  values  and  comparing  them  with  values 
from  formula  2,  the  approximate  values  using  this  formula  have  been  found  to 
be  very  accurate  for  typical  earth  values  and  frequencies  in  the  HF  range. 

Some  typical  values  for  the  line  wavelength  and  attenuation  using  this  second 
formula  for  are  shown  in  table  1.  Also  shown  for  air  insulation  at 
4  and  16  MHz  for  ground  parameters  of  permittivity  16  and  conductivity  .01  S/M. 

IABLE  1 


4  MHz 

16  MHz 

Cladding 

Aline 

Abeat 

loss 

Aline 

Abeat 

loss 

thickness 

metres 

metres 

dB/M 

metres 

metres 

dB/M 

none 

13.28 

16.14 

2.9 

4.44 

5.83 

3.81 

•  5m 

17.41 

22.68 

.79 

5.76 

8.32 

1.93 

1  mm 

21.41 

29.98 

.46 

6.66 

10.3 

1.45 

2  ran 

26.57 

41.17 

.28 

7.85 

13.5 

1.05 

5  ran 

34.28 

63.17 

.16 

9.70 

20.11 

.719 

10  ran 

40.17 

86.58 

.12 

11.15 

27.56 

.574 

25  ran 

47.54 

130.0 

.09 

13.00 

42.51 

.455 

infinity 

74.98 

inf. 

.00 

18.74 

inf. 

.000 

Since  the  size  of  the  insulation  is  a  small  fraction  of  a  wavelength  it 
would  be  expected  that  higher  order  modes  would  be  rapidly  attenuated. 

Consider  an  antenna  such  as  shown  in  figure  1  driven  at  one  end  of  the 
insulation  against  an  earth  spike.  The  first  part  of  the  line  will  set  up  a 
lateral  wave  travelling  along  the  surface  with  a  wavenumber  K  (of  free  space) 
but  since  the  wave  is  travelling  along  the  line  with  wavenumber  Kj,  the  relative 
phase  of  the  line  and  the  lateral  wave  goes  up  and  down  with  a  wavelength 
beat  given  in  the  table  where 


l/X, 


beat 


1/Xline“1/Xair 


This  effect  is  clearly  shown  in  figures  2  and  3.  These  figures  show  the 
relative  received  power  at  the  surface  for  the  stated  length  of  buried 
insulated  antenna.  It  is  assumed  a  unit  wave  has  been  launched  down  the 
antenna  and  the  receiving  point  is  on  the  surface  in  a  vertical  plane 
containing  the  antenna  in  a  direction  away  from  the  antenna  the  same  as  the 
direction  of  the  travelling  wave. 
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Values  plotted  in  figures  2  and  3  are  calculated  from 


<|/*X  -  K.irldl|  !/<|/el(KL  -  Kair>dl  > 

Where  K.  is  the  wavenumber  of  the  line,  K  .  the  wavenumber  in  air. 
l.  air 

This  is  the  signal  relative  to  a  1  metre  length. 

The  thicker  the  cladding  the  larger  the  eventual  signal  will  be,  since 
the  longer  the  length  of  line  before  phase  reversal  occurs.  After  the 
maximum  is  reached  the  signal  falls  away  again.  The  effect  of  the  attenuation 
in  the  line  is  to  damp  the  oscillation  shown  and  the  signal  eventually  fallB 
to  a  constant  value  unaffected  by  any  change  in  length.  In  the  antenna  is 
short,  say  for  example  10  metres,  and  since  all  thicknesses  of  cladding  have 
the  same  efficiency  curve  at  first,  it  would  be  pointless  to  increase  the 
insulation  beyond  a  certain  thickness  indicated  by  these  curves.  The  effect  of 
a  plastic  insulator  instead  of  an  air  one  is  seen  to  decrease  the  beat  wave 
length  and  so  decrease  the  size  of  the  eventual  signal  maximum. 

The  front  to  back  ratio  can  also  be  calculated  in  a  similar  simple  way 
and  is  shown  in  figures  4  and  5.  Short  lengths  do  not  shown  any  back  to  front 
ratio. 


The  backward  signal  behaves  in  a  manner  similar  to  the  forward  signal  but 
since  the  phase  of  the  wave  cn  the  line  and  lateral  wave  are  in  opposition 
the  beat  wavelength  is  much  more  rapid  and  is  given  by:- 


l/Xbe.t  * 


Again  due  to  attenuation  there  is  not  complete  cancellation  and  the 
backward  signal  falls  to  a  small  value  which  is  not  zero.  The  front  to  back 
ratio  therefore  rises  cusp-like  to  a  maximum  at  the  first  minimum  of  the 
backward  signal  given  by  the  above  formula. 

The  effect  of  depth  of  burial  can  be  assessed  by  noting  that  for  a 
reasonable  earth  dielectric  constant  the  angle  of  travel  to  the  surface  is 
nearly  vertical  and  these  rays  suffer  an  attenuation  given  in  table  1  of 
2.9  dB/m  at  4  MHz  and  3.88  dB/M  at  16  MHz  for  the  ground  parameters  given. 

3  MOMENT  METHOD  TREATMENT 


In  the  previous  transmission  line  treatment  there  are  a  number  of  effects 
that  are  ignored  ie  the  effect  of  the  earth's  surface  on  the  transmission  line 
and  the  practical  excitation  of  the  line  itself  which  may  be  uncertain,  as 
there  is  no  definite  outer  conductor  to  connect  to.  To  assess  the  accuracy  of 
the  transmission  line  treatment  a  moment  method  calculation  was  carried  out. 

As  is  usual  with  the  moment  method  [  4]  the  central  conductor  was  divided  into  a 
number  of  elements  along  its  length.  Each  one  is  assigned  an  unknown  current. 
The  electric  field  at  one  element  due  to  the  current  in  another  is  calculated 
using  the  routines  in  reference  [  5] .  In  addition  the  effect  of  the  dielectric 
cladding  was  included  by  a  method  due  to  Richmond  [  6] ,  described  more  fully  in 
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an  appendix.  The  boundary  condition  of  zero  tangential  electric  field  at  the 
conductor  surface  was  imposed  at  the  centre  of  each  element.  Solving  the 
resultant  set  of  simultaneous  equations  enables  the  current  on  the  conducting 
wire  to  be  calculated.  Knowing  this  current  the  field  on  the  surface  of  the 
earth  at  a  distance  can  be  calculated  using  the  methods  described  in 
reference  {  7] .  This  essentially  consists  of  adding  the  attenuation  of  a  wave 
travelling  to  the  surface  to  ther  attenuation  of  the  lateral  wave  travelling 
along  the  surface. 

The  current  in  the  central  conductor  is  shown  in  figure  6,  the  bottom 
part  of  the  figure  shows  schematically  the  arrangement  of  conductors  and 
insulation.  There  is  a  central  insulated  part  of  20  metres  in  length  and  two 
bare  sections  of  10  metres  each  end  which  serve  as  earth  spikes.  The  antenna 
is  driven  by  a  unit  voltage  generator  at  the  junction  of  the  left  hand  bare 
section  with  the  insulated  part  ie  at  the  10  metre  point  on  the  length  scale. 

The  two  curves  above  this  show  the  relative  amplitude  and  phase  respectively 
of  the  current  in  the  conductor.  To  the  left  of  the  driving  point  the 
current  decays  exponentially  fairly  rapidly  with  a  phase  that  changes  linearly 
with  distance.  This  clearly  demonstrates  an  attenuated  travelling  wave  down  the 
earth  spike.  The  complex  wavenumber  for  this  wave,  if  the  wire  is  thin,  will 
be  the  complex  wavenumber  for  waves  travelling  in  the  conducting  dielectric 
earth.  To  the  right  of  the  driving  point  another  wave  travels  with  less 
attenuation  and  a  less  rapid  phase  change  with  distance  down  the  insulated 
part  of  the  wire,  thus  exhibiting  transmission  line  behaviour.  Note  the 
ripples  in  amplitude  at  the  end  of  the  insulated  part  indicating  a  reflection 
at  the  end  and  a  reflected  wave  travelling  back  down  the  line.  In  this  case 
the  amplitude  reflection  factor  deduced  from  the  ripple  amplitude  is  about  15%. 

From  the  plots  of  phase  and  relative  amplitude  of  the  current  the  line 
wavelength  and  attenuation  given  by  the  moment  method  can  be  compared  with  the 
transmission  line  calculation. 


TABLE  2 


Moment  Method 

Transmission 

line 

FREQ 

Wavelength 

Loss 

Wavelength 

Loss 

MHz 

Metres 

dB/M 

Metres 

dB/M 

4 

32. 

- 

34.28 

.164 

11 

13.7 

.61 

13.6 

.60 

16 

9.8 

.81 

9.7 

.719 

The  two  methods  agree  quite  well. 

Also  printed  at  the  bottom  of  figure  6  is  the  field  in  the  forward 
direction  (along  the  antenna  from  left  to  right)  on  the  surface  of  the  earth  v 
at  a  range  of  1000  metres.  The  front  to  back  ratio  and  impedance  at  the 
driving  point  are  also  given. 

A  15%  reflection  at  the  insulation  end  could  degrade  an  infinite  front  to 
back  ratio  to  16.5  dB.  The  reflection  at  the  end  is  therefore  only  really 
important  for  high  front  to  back  ratios,  which  as  can  be  seen  from 
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figures  4  and  5  occurs  only  for  narrow  ranges  of  antenna  length.  This 
reflection  at  the  end  of  the  insulation  can  be  reduced  by  matching  the 
insulated  transmission  line  part  to  the  bare  part.  A  series  of  resistances 
were  found  to  be  ineffective  as  they  tended  to  be  shunted  by  the  conducting 
earth.  Another  method  shown  in  figure  7  is  to  taper  the  insulation  down  to 
zero.  As  can  be  seen  from  the  plots  of  the  amplitude  and  phase  of  the  current 
the  background  reflected  wave  has  been  reduced.  Due  to  the  relatively  long 
length  of  the  antenna  the  front  to  back  ratio  is  only  6.6  dB  in  figure  6  and 
reducing  the  backward  reflection  would  have  little  effect.  Increasing  the 
effective  length  of  the  antenna,  a  much  more  critical  factor  has  in  fact 
reduced  the  front  to  back  ratio. 

The  left  hand  bare  section  or  earth  spike  of  figure  7  shows  a  rapid 
exponential  decay  to  a  small  value,  the  shorter  earth  spike  shown  in  figure  8 
is  adequate.  The  bare  section  at  the  end  of  the  tapered  section  is  still 
needed  as  is  evident  from  figure  9. 

Since  the  moment  method  calculations  exhibit  transmission  line  behaviour 
it  would  be  expected  that  such  calculations  agree  with  those  of  figures  2-5. 

To  confirm  this  a  set  of  different  lengths  were  analysed  at  4  MHz,  one  of  the 
set  is  shown  in  figure  10.  Here  it  is  seen  that  the  taper  is  only  partially 
effective  in  reducing  the  back  reflection,  but  as  previously  explained  this  is 
only  important  for  high  back  to  front  ratios. 

Table  3  gives  the  front  to  back  ratio  as  a  function  of  the  total 
insulator  length  including  the  tapered  section. 


TABLE  3 


Length 

F/B  ratio 

Gain  relative 

metres 

dB 

to  10  m  in  dB 

10  m 

1.4 

0 

14  m 

3.3 

3.64 

20  m 

13.1 

4.24 

24  m 

16.4 

3.52 

These  results  compare  quite  well  with  figures  2  and  4  the  transmission 
line  calculations. 


The  4  metre  tapered  section  is  about  .3  transmission  line  wavelengths  at 
11  MHz  and  appears  adequate,  whereas  at  4  MHz  it  has  fallen  to  .12  wavelengths 
and  is  clearly  too  short.  This  suggests  a  quarter  wavelength  is  about  the 
minimum  length. 

Results  at  2  MHz  require  a  longer  taper  and  a  set  of  results  was 
calculated  which  illustrates  the  effect  of  changes  in  taper  length  and  also 
the  effect  of  earth  spikes.  The  transmission  line  wavelength  is  about  64  metres 
for  the  2  MHz  case,  so  that  the  taper  should  be  at  least  16  metres.  Table  4 
gives  front  to  back  ratios  as  a  function  of  length  and  taper  length:- 
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Figure  11  shows  an  antenna  with  44  ro  of  insulation  at  2  MHz,  there  is  an 
appreciable  reflection  but  despite  this  the  front  to  back  ratio  is  still 
13.1  dB.  The  reflection  from  the  end  is  reduced  by  the  10  m  taper  shown  in 
figure  12.  Figure  13  is  an  amplitude  and  phase  diagram  for  this  antenna. 
Starling  at  the  origin  a  small  vector  is  added  corresponding  to  each  current 
element  along  the  line,  one  diagram  for  the  forward  and  one  for  the  backward 
direction.  The  section  from  the  origin  to  the  first  slanting  arrow  is  the 
signal  from  the  first  bare  part  or  earth  spike.  The  section  to  the  next  arrow 
is  the  signal  from  the  5  mm  thick  insulated  part,  the  section  between  the 
second  and  third  slanting  arrows  is  the  taper  and  finally  the  signal  from  the 
second  (right  hand  side)  earth  spike  which  is  negligeable.  Although  the 
backward  signal  diagram  is  far  from  the  ideal  circle  the  resultant  backward 
signal  is  not  large. 

In  figure  14  the  taper  is  longer  so  reducing  the  backward  signal  in  the 
transmission  line  but  due  to  the  signal  from  the  earth  spike  at  the  beginning 
of  the  insulated  part  the  resultant  backward  signal  is  small  and  the  front  to 
back  ratio  large. 

Figure  15  shows  a  longer  taper  which  has  actually  reduced  the  front  to  back 
ratio  as  the  amplitude  and  phase  diagram  for  the  backward  signal  has  wound  up 
too  tightly  leaving  the  signal  from  the  first  earth  spike  remaining.  Running 
this  earth  spike  at  right  angles  to  the  antenna  could  reduce  this  remaining 
signal.  Reducing  the  backward  signal  therefore,  may  not  always  be  a  good  thing. 

Table  4  summarises  the  2  MHz  results.  The  variation  of  the  front  to  back 
ratio  for  antennas  with  14  metre  tapers  ie  the  24,  34  and  44  metre  lengths 
agree  quite  well  with  the  transmission  line  calculations  at  the  top  left  of 
figure  4. 


TABLE  4 


Total  Insulation 
length  metres 

Taper  length 
metres 

F/B 

dB 

Field 

1  rnW 

at  1000  m  for 
drive  power 

24 

14 

3.0 

5.6 

34 

14 

8.2 

6.9 

44 

0 

13.1 

6.0 

(Fig  11) 

44 

10 

19.7 

6.1 

(Fig  12,  13) 

44 

14 

25.1 

6.1 

(Fig  14) 

44 

20 

15.8 

6.0 

(Fig  15) 

4  CONCLUSIONS 

Both  transmission  line  methods  and  moment  methods  can  be  used  to  analyse 
buried  antennas.  Factors  relevant  to  buried  antenna  design  are:- 

4 . 1  LENGTH 

Reference  to  figures  1  to  4  indicates  what  front  to  back  ratio  and 
relative  gain  can  be  expected  for  various  lengths  and  frequencies. 
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4.2  INSULATION  THICKNESS 


Given  a  value  for  the  length  it  can  be  seen  from  figures  1  to  4 
there  is  a  thickness  above  which  it  is  not  worth  increasing  the 
insulation  thickness. 

4.3  DEPTH  OF  EURIAL 

The  loss  is  given  in  table  1  and  is  around  a  few  dB  per  metre  of 
burial.  This  would  not  be  significant  for  shallow  burial. 

4.4  LENGTH  OF  EARTH  SPIKES 

Again  table  1  gives  the  attenuation  for  a  bare  wire  which  is  quite 
high.  The  length  should  be  chosen  to  obtain  a  reasonable  attenuation 
say  80%.  A  length  longer  than  this  is  not  worthwhie. 

4.5  TAPERING  OF  THE  INSULATION 

The  signal  reflected  off  the  end  of  the  effective  transmission  line 
part  of  the  antenna  can  be  reduced  by  tapering  the  insulation  thickness  to 
zero . 

4 . 6  PLASTIC  INSULATION 

This  is  appreciably  worse  than  the  air  insulated  case  from  the  gain 
and  front  to  back  ratio  point  of  view,  as  can  be  seen  from  figures  2  to 
5. 
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APPENDIX 


THE  RICHMOND  DIELECTRIC  CLADDINC  APPROXIMATION 

This  1 s  described  in  reference  6.  The  wire,  cladding  and  surrounding 
earth  are  shown  in  figure  Al.  The  first  step  is  to  rearrange  Maxwells 
equation  for  inside  the  cladding  region  of  dielectric  constant 

Curl  H  -  J  +  D 

=  J  ♦  "  ex  +  c^E 

=  J  +  CiE  +  (e2  -  Cl)E 

it  the  cladding  can  be  replaced  by  an  equivalent  current 

J  .  =  jaj(e -  -  r,  )E 

equiv  J  2  1 

assuming  a  time  dependance  of  e^Wt 

near  a  wire  the  field  is  perpendicular  to  the  wire.  The  charge  on  a  length 

...  .  3 1  dl 

•j  ■.  or  wire  is  -  —  x  — 

a  5  jx 

the  racial  field  in  the  direction  S  is  given  by 


'RADIAL 


e2ju;.  2-rrS 


and  so  J 


equiv 


—  •  (-aA  i 

\  3  5.  /  S 


( t;2  ~  CI) 


-  Al'  —  where  A  =  - 

S  2ne 


I'-^l 


referring  to  Figure  A2 

the  current  in  the  sector  of  angle 

=  Al'dfa? 

for  a  disc  $  -  2tt 
Total  Current  =  27rAI'd£ 


knowinb  this  current  the  field  at  any  point  due  to  it.  may  be  calculated,  note 
it  is  independent  of  S  the  radius. 
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